Pasta phases in neutron star studied with extended relativistic mean field models 
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To explain several properties of finite nuclei, infinite matter, and neutron stars in a unified way 
within the relativistic mean field models, it is important to extend them either with higher order 
couplings or with density-dependent couplings. These extensions are known to have strong impact 
in the high-density regime. Here we explore their role on the equation of state at densities lower 
than the saturation density of finite nuclei which govern the phase transitions associated with pasta 
structures in the crust of neutron stars. 
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I. INTRODUCTION 

The outer layer of the neutron star (NS), with density 
less than the nuclear saturation density, represents dif- 
ferent challenges and observational opportunities such as 
thermal evolution, x-ray burst, glitches, and the very im- 
portant core-crust transition region At this density, 
nucleons are correlated at short distances by attractive 
strong interactions- they are anticorrelated at large dis- 
tances because of the Coulomb repulsion. Competition 
among short- and long-range interactions leads to the de- 
velopment of complex and exotic nuclear shapes, which 
can be oversimplified @ to spheres, bubbles, rods, slabs, 
and tubes. These are generally known as pasta phases. 
Similar geometries can also be obtained for uncharged 
systems owing to finite size effects Q. In NS matter, ge- 
ometries such as slab with holes cross- rods ||, jungle 
gym, curled spaghetti @, etc., could be possible but are 
neglected in the present Brief Report. The pasta phases 
in the crust eventually dissolve into uniform matter at 
a certain larger density close to the saturation density. 
Such a critical density at which the mantle-crust transi- 
tion in NSs happens can be related to the neutron skin 
thickness in heavy nuclei Q ■ Several other links between 
the NSs and finite nuclei are established (8l-[To| and hence 
formulating unified models which explain both finite nu- 
clei and NSs become important [uj. 

The relativistic mean field (RMF) models are success- 
ful in explaining several properties of finite nuclei. The 
equation of state (EoS) for symmetric matter and NS 
matter obtained from many such standard interactions 
do not satisfy the experimental/observational constraints 
@, QI3, [HI i mainly because of their inconsistent behavior 
at high density. In such cases, it is useful to consider ex- 
tensions in the RMF models carried out by formulating 
an effective hadron field theory either with additional 
couplings @, [HI- or with density-dependent meson nu- 



cleoli couplings [llj. These extensions do not have a 



strong impact on the EoS near saturation densities (be- 
cause the parameters are fitted with data for finite nuclei) 
but they strongly modify the EoS at densities higher than 
the saturation density. In the present Brief Report, we 
analyze the impact of the extension in the RMF models 
on the EoS at lower densities which govern the occurrence 



of pasta phases. 

In case of higher order couplings one may relate this 
analysis to a trivial fact that if we fit few data points 
with higher order polynomials, the extrapolated results 
in both of the sides away from the data strongly depend 
on the order of the polynomial chosen. In the case of 
EoS, the energy density decreases with the decreasing 
density to reach zero at zero density, and so are the cor- 
responding fields/interactions. Hence, the effect of higher 
order couplings are not expected to be drastic but could 
be appreciable at finite subsaturation densities. Pasta 
structures studied with higher order couplings were first 
reported by us [l2| and a more comprehensive analysis 
of different nature has recently been published [l3j| . 



II. FORMALISM 

In this section we briefly discuss the formalism for cal- 
culating the energy of nuclear matter in homogenous and 
pasta phases. For the homogenous phase, we employ 
RMF models with parameters of different classes namely,, 

1. standard nonlinear interaction (NL3 [Hj|); 

2. density dependent interaction (DDME2 [U); 

3. interactions with a few additional hig her order cou- 
plings (FSUGold and IUFSU to constrain 
selected observables; and 

4. Effective field theory motivated interactions with 
several higher order couplings (Gl and G2 fl6j). 

These effective interactions are adjusted to reproduce 
various properties of selected finite nuclei and tested very 
well by explaining experimentally observed features in 
a variety of nuclei. The details of the Lagrangian and 
the coupling constants can be found in the correspond- 
ing references given above. The energy density for the 
homogenous nuclear matter is given as 
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FIG. 1: (Color online) Energy per particle calculated with 
different parameter sets, for proton fractions 0.1, 0.3, and 
0.5. The numbers adjacent to the name of the parameters 
represent the corresponding neutron skin thickness (in fm) 
obtained for 208 Pb. 
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where k fp = (S^xp) 1 / 3 , k fn = (37r 2 (l - x)p) 1 ^ 3 , p is the 
total baryon density, and x defines the proton fraction. 
In the above expressions a, Vq, and Rq denote the scalar, 
vector, and isovector meson mean fields, m a , m u , and 
trip are the corresponding meson masses, and m„ is the 
nucleon mass. The symbols g a , g^, g p , k 3 , k 4 , rji, r] 2 , r\ p , 
A„ and Co denote the various coupling constants. 

For the pasta phase, we employ the Wigner-Seitz ap- 
proximation to calculate the Coulomb and surface en- 
ergies for the different pasta structures. Wigner-Seitz 
cells arc approximated by (i) a sphere of radius r w in the 
case of spheres and bubbles, (ii) a cylinder of radius r w 
in the case of rods and tubes, and (hi) a cell bounded 
by planes (with r w being defined as the half distance be- 
tween plane boundaries) in the case of slabs P, Q ■ 

Then total energy density (EoS) of the pasta phase Q 

is 



(2) 



where e e , e c , and e s represent the electron, surface, and 
Coulomb energy densities, respectively, and it represents 
the ratio of the volume of the cluster and the Wigner- 
Seitz cell. In other words, u represents the filling factor 
with u = n/n', where n and n' represent the average and 
dense phase density, 

e s = uad/r and e c = 2u{e 2 TT(f) 1/3 {(jdxn' /4) 2 / 3 , (3) 

where u = u for droplets and u = 1 — u for bubbles, a 
is the surface tension coefficient calculated by a Thomas- 



300 - 



> 

CD 



60 - 



— i 1 1 1 1 1 1 1 1 1- 

(a) x=0.1 




■ — I — ■ — I — I — I — I — I — ■ — I 



(c) 



FSUGold 

IUFSU 

DDME2 




0.00 0.02 0.04 0.06 0.08 0.10 0.12 

Density (fm" 3 ) 

FIG. 2: (Color online) Density dependence of the scalar (ffcer), 
vector (gtaVo), and isovector {g P Ro) fields in the NS matter 
calculated with different parameter sets and at the proton 
fraction x = 0.1 . 



Fermi method j!7| . d — 1,2,3 represents the dimension 
of the system, and 



1 



(d + 2) 



(4) 



In the case of d = 2, the above expression can be rewrit- 
ten as 



(u — 1) — log u 



(5) 



The radius of the droplet (rod, slab) and that of the 
Wigner-Seitz cell are, respectively, given by 



r/u 



l/d 



(6) 



III. RESULTS AND DISCUSSION 



We first examine the variation of energy at subsatu- 
ration densities relevant for pasta phase. Figure [T] rep- 
resents the energy per particle for nuclear matter as a 
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FIG. 3: (Color online) 
Sphere Transition between 
p oc j different pasta config- 
Slab urations, for proton 
Tube fractions 0.1, 0.3, and 
Bubble 0.5, calculated with 
Uniform different parameter sets. 
The different colors 
represent sphere, rod, 
slab, tube, bubble and 
uniform phases as given 
in the legend. 



function of density calculated with different parameter 
sets and proton fractions. We can see the general trend 
that with increasing proton fraction the matter is becom- 
ing more and more bound. For symmetric nuclear matter 
(proton fraction x — 0.5) the binding is maximum and 
we do not sec much difference in the results from different 
parameters. The differences start to appear at x = 0.3 
and with x = 0.1, we can identify three clear groups 
of parameter sets: (i) having a shallow well in the en- 
ergy (G1,G2, NL3), (ii) saturated energy at low density 
(FSUGold, DDME2), and (iii) with hump in the energy 
at low density (IUFSU). To understand the above group- 
ing and related features in terms of the different fields 
and their couplings, we look into the density dependence 
of the a. ui, and p fields, which are shown in Fig. [2] We 
notice that at subsaturation densities, the a field is the 
most dominant. The higher order terms in the a field 
(self-coupling) are taken care of in all the chosen param- 
eters upto the fourth order. DDME2 is an exception 
with linear a field but with explicitly density-dependent 
couplings which decrease drastically with density for the 
range shown in Fig. [2] The u> field follows the same pat- 
tern of the a field and hence we may not be able to see 
the explicit role of higher order couplings involving the 
u> field alone or those involving both a and uj. 

We see quite interesting features in the p field which 
is almost linear in Fig. M.c) for NL3, Gl and G2. The 
difference in energy owing to the change in p field could 
be expected through the second, fifth and sixth terms in 
Eq. (p}. Among these terms the second and fifth terms 
can contribute to the nonlinear behavior of p field (while 
a and uj fields are almost linear and the proton frac- 
tion is constant). Between these two terms, the second 
one (with the coupling constant A„) dominates but this 
term appears exclusively in the parameters FSUGold and 
IUFSU. Thus the higher order coupling A„, representing 
the second order uj-p interaction, plays a major role and 
we can identify that the hump in the energy is propor- 
tional to A„. This parameter plays an important role at 
higher densities also, as evident from its influence on NS 
mass and radius 1181 . 



The parameter set DDME2, even without the higher 
order coupling, yields a p field similar to that of FSUG- 
old because the coupling constant g p varies nonlinearly 
with density [ll[ and mimics the role of A„ . It has to be 
noted that the neutron skin thickness has been consid- 
ered to be an ingredient in the fitting procedure of the 
three parameters DDME2, FSUGold and IUFSU. The 
values of neutron skin thickness obtained for 208 Pb are 
also shown in Fig. [TJ along with the legends. We can see 
that the ordering of energy at x = 0.1 for different pa- 
rameters is same as that of the neutron skin thicknesses. 
As the neutron skin thickness is strongly correlated to 
the density-dependence of symmetry energy, we can as- 
sociate the differences in energies to the differences in 
the symmetry energies. Hence it is more appropriate to 
link the low-density behavior of EoS (especially for neu- 
tron rich case) to the symmetry energy as discussed in 
Ref . [TtJ ■ The surface tension coefficients [l?} also follow 
the same ordering of the energy and neutron skin thick- 
ness. In Fig. 1, we can see that the different parameters 
converge at a density lower than the saturation density 
(« 0.11 fin" 3 for x = 0.1). This is attributable to the 
fact that the parameters arc fitted to reproduce the prop- 
erties of finite nuclei whose average density is lower than 
the saturation density (l9j . 

We identify that the p field is the major source of 
the difference in the EoS from different parameters. At 
higher proton fractions the a and uj fields remain almost 
the same and the strength of the p field decreases. Thus, 
the results from different parameters converge at higher 
proton fraction. 

The various phase transitions happening in the NS 
comprising the considered pasta structures, obtained 
from our calculations with different proton fractions and 
parameter sets, are represented in Fig. [3] It has to be 
noted that there could be more geometries of pasta struc- 
tures 4- 6] and hence Fig. [3] is not a complete phase di- 
agram. With all the parameter sets, we find almost the 
same trend in the sequence of pasta structures. At low 
density, sphere phase is the most stable configuration, fol- 
lowed by rod, slab, tube, and bubble phases. In most of 
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FIG. 4: (Color online) Energy per 
particle for the uniform matter and 
pasta-like matter obtained with the 
IUFSU parameter set for proton 
fraction (a) 0.1, (b) 0.3, and (c) 0.5. 
The different colors represent the 
different pasta structures as men- 
tioned in the inset. The energy con- 
tribution from the electrons is also 
included here (unlike Fig. 1). 
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the cases, a homogeneous matter phase occurs at a den- 
sity ~ 0.11 fm~ 3 , which, incidentally, is the point where 
all the EoS converge [lj| . The difference between results 
from different parameter sets could be clearly linked to 
the features discussed above in the case of EoS [Fig. 1]. 
Similar to the case of energy, at higher proton fractions 
the different parameter sets yield similar results and we 
see interesting features at x = 0.1. The sustained ex- 
istence of uniform phase at lower densities in the cases 
of FSUGold, IUFSU and DDME2 can be linked to the 
corresponding EoS shown in Fig. [TJ which are larger and 
raising. The energy difference between uniform phase 
and pasta phase is very small at x = 0.1, and such a 
delicate balance allows the differences in EoS to show up 
significantly by altering the transition densities from one 
phase to other. 

In Fig 4, we compare the energies calculated for the 
uniform phase and the pasta phases, in case of IUFSU 
for which the difference with other parameters is notable 
in Fig. 3. The crossing between these energies determine 
the corresponding phase transitions. The main differ- 
ence between the uniform phase energy in Fig. [1] and 
Fig. H is that in Fig. H we include the contribution from 
electron energy also. In general, the pasta phase energy 
[Eq.ljJJ] is dominated by the corresponding uniform phase 
energy. At x = 0.1, the contribution from the Coulomb 



term is lesser and hence the difference between energies 
of pasta phases and uniform phase is smaller, as evident 
from Fig. 4. At larger proton fractions the Coulomb term 
widens the energy gap and the difference between dif- 
ferent EoSs starts to diminish, leading to similar phase 
transitions. 



IV. CONCLUSIONS 

With lower proton fractions, the higher order couplings 
can play a role in modifying the EoS at subsaturation 
densities relevant for the pasta structures in neutrons 
stars. The a field dominates the EoS at subsaturation 
density and its higher order couplings are taken care of 
all the realistic models. The fourth-order coupling be- 
tween uj and p fields (A„), plays a crucial role at lower 
proton fraction and a similar role is played by the density- 
dependent p field coupling constant. At higher proton 
fractions, the p field, which is the source of the difference 
between different parameter sets, gets weaker and hence 
all the results converge. We observe that it is more ap- 
propriate to link the differences in EoS (and hence the 
occurrence of pasta structures) to the symmetry energy 
(and hence the neutron skin thickness). 
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